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ABSTRACT 



Context. The investigation of the dominant formation mechanism of early-type galaxies. 

Aims. In this work we derive clues to the formation scenario and ionisation source of the lenticular galaxy NGC 3607 
by means of metallicity gradients, stellar population and emission lines properties. 

Methods. We work with long-slit spectroscopy from which we (i) study the radial distribution of the equivalent widths 
of conspicuous metallic absorption features, (ii) infer on the star-formation history (with a stellar population synthesis 
algorithm), and (iii) investigate the ionisation source responsible for a few strong emission lines. 

Results. Negative radial gradients are observed for most of the absorption features of NGC 3607. Compared to the 
external parts, the central region has a deficiency of alpha elements and higher metallicity, which implies different 
star-formation histories in both regions. At least three star formation episodes are detected, with ages within 1-13 Gyr. 
The dynamical mass and the Mg-2 gradient slope are consistent with mergers being important contributors to the 
formation mechanism of NGC 3607. Emission-line ratios indicate the presence of a LINER at the centre of NGC 3607. 
Contribution of hot, old stars to the gas ionisation outside the central region is detected. 

Conclusions. Evidence drawn from this work suggest small mergers as important contributors to the formation of 
NGC 3607, a scenario consistent with the star-formation episodes. 

Key words. Galaxies: elliptical and lenticular, cD; Galaxies: individual: NGC 3607; Galaxies: fundamental parameters 



1. Introduction 

Significant efforts have been undertaken in the last few 
years to understand the complex mechanisms associated 
with the formation and evolution of the early-type galaxies. 
Since the early 90's it was established that ellipticals cannot 
be described as being part of a one-parameter family, and 
that many show unambiguous signat ures of interaction wit h 
companions and the environment (jvan den Berghl fl990). 
Such facts suggest that these galaxies may have experi- 
enced different star formation histories, with the present- 
day stellar populations differing in metallicity and/or age. 
In this sense, most early-type galaxies (or their nuclei) 
are not characterised by a single-aged stellar population. 
Instead, they are better described by a composition of single 
stellar populations (e.g. TWorthev. Faber fc Gonzaleslf l992: 
Arimoto fc Yoshiil Il987t iRickes et alJ 12004 iRickes et al I 



2008h . 



Observationally, clues to the star formation history of a 
galaxy can be found in the radial dependence of the stel- 
lar [q /Fel ratio, or more specifically the M g/Fe abundance 
ratio (jPipino. Matteucci & Chiappinil[200l ). The Mg2 line- 
strength distribution in early-type galaxies can vary consid- 
erably, with radial profiles ranging from essentially feature- 
less to rather structured. In the latter case, changes of slope, 
possibly associated with kinematically decoupled cores, or 
anomalies in the stel lar population, hav e been observed in 
some ellipticals (e.g. ICarollo et aLlll993f ). 



In this context, the investigation of the spatial distri- 
bution of the stellar population and metallicity in samples 
of early-type galaxies may shed some light on the forma- 
tion and evolution issues. With the above goals in mind, we 
have been conducting a systematic stud y of objects rich in 
interstellar med ium, such as NGC 5044 (jRickes et al.| [2004) 
and NGC 6868 (IRickes et al.ll2008h . 

This work focus on the lenticular galaxy NGC 3607, 
which is the central and brightest member of a group that 
contains the galaxies NGC 3605 and NGC 3608 as well. 
Fig. [TJ shows a 15' x 15' V band image of NGC 3607 and 
companions, extracted from the Canadian Astronomy Data 
Centre (CADC0 

The radial velocity of NGC 3607 is - 960 kms" 1 
([Ferguson &: Sandagd Il990l ) , which puts it at a distance 
of <~ 12.8Mpc (with H = 75kms~ 1 Mpc" 1 ). At such a 
distance, 1" corresponds to a bout 0.062 kpc. NGC 36 07 is 
a LINER with X-ray emission (jTerashima et al.1 12002) , and 
discrete sources in the host galaxy are the probable origin of 
the extended hard X-ray emission . The hot gas in th i s ellip - 
tical galaxy has been analysed bv lMatsushita et al.l ((2000). 
NGC 3607 is another case where the dust and ionised gas 
distribute as a small asymmetric disk. The dust absorp- 
tion is stronger in the northeast direction, while the ionised 
gas emission is m ore important towards the southwest 
([Ferrari et aTlll999D . 

The goal of this paper is to investigate the stellar pop- 
ulation, metallicity distribution and the presence of ionised 
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Table 1. General data on NGC3607 




Fig. 1. V band image covering the field 15' x 15' centred on 
NGC 3607. Image orientation: North to the top and East 
to the left. 

gas in NGC 3607, which are important parameters to un- 
derstand the formation and evolution of this galaxy. 

This paper is organised as follows: Sect. [2] presents the 
observations and data reduction; Sect. [3] deals with the 
equivalent width analysis of the absorption lines and their 
radial dependence; Sect. 4 describes the stellar population 
synthesis, whose results are discussed in Sect. 5; Sect. 6 
presents an analysis of the emission lines and the nature of 
the ionisation source. Finally, in Sect. 7 a general discussion 
of the results and conclusions will be drawn. 



2. Observations and spectra extraction 

In the present paper we study the star formation his- 
tory, metallicity distribution, and emission-gas properties 
of NGC 3607. Basically, we work with long-slit spectroscopy 
and a star formation synthesis algorithm, following simi- 
lar methods as tho se applied to NGC 6868 and NGC 5903 
(Ric kes et al.l l2008h . General parameters of the bright SO 
galaxy NGC 3607 are given in Table [Q 

The observations were obtained with the ESO 3.6m tele- 
scope at La Silla, Chile, equipped with EFOSC1. The spec- 
tra cover the range AA5100 — 6800 A with a 3.6A/pixel res- 
olution. The spatial scale of the observational configuration 
is 0.6" pixel -1 ; the slit length was 3.1' while its width was 
fixed at 1.5", corresponding approximately to the average 
seeing; two 2-dimensional spectra were obtained with the 
same exposition time (approximately 45 minutes) and po- 
sition angle PA = 301°. In all cases, the extraction area 
was fixed at 4.57" 2 , which corresponds to « 0.02 kpc 2 . 

The standard star HR 7333 was observed for flux and 
velocity calibrations. We point out that the average FWHM 
of the absorption features in the spectrum of HR7333 is 
between 8 and 10 A, similar to that of the Lick/IDS system. 

The spectra were corrected for the radial velocity 
V r = 960kms _1 and the foreground Galactic extinction 



Parameter 



Value 



a(J2000) 
5(J2000) 

Mb 
B 

E(B - V) 

a+[N II] 

Heliocentric radial velocity 

Distance 

Redshift 

Diameters 

Effective radius 1 



iri6 m 54.6 s 
+ 18°03'07" 
-21.92 
10.8 
0.070 

1.28 x 10 41 ergs" 1 
960 ± 20 km s" 1 
12.8 Mpc 

0.003202 ± 0.000067 

4.9' x 2.5' 

43.4" 



Table Notes. Data obtained at the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under 
contract with the Natio nal Aeronauti c s an d Space 
Administration, (f) R e from I Annibali et all (|2007l ). 

Ay = 0.035 (Table QJ. The corrected spectra, normalised 
at A5870A, are shown in Fig. [2j separately for the northeast 
and southwest directions on the galaxy. 

3. Radial properties of absorption features 

Important information on the chemical abundance and stel- 
lar population of a galaxy can be obtained from the analysis 
of the equivalent widths of a set of selected absorption lines 
and c ontinuum fluxes (e.g. iRickes et alJ 120081 : [Rickes et al.l 
|2004| . and references therein). Here we work with the con- 
tinuum tracing and spectral windo ws for the absorp tion 
features defined by the Lick system (Fa ber et al.lll985l ). 

Because of the somewhat restricted coverage of our 
NGC 3607 spectra (Sect. [2]), the available absorption fea- 
tures are Mg2, Fei A52 7o, Fei A5 3 35 , Fei A54 o 6 , Fei A5709 , 
Fel A578 2, Nai A589 5, and TiOn A6 237- Continuum points 
within the same spectral region are A5300, A5546, A5650, 

A58OO, A5870, A6173, A6620, and Ag640- 

Some early-type galaxies have been shown to present 
systematic spatial variations on the line strength of 
the indices Mg2 and Fei, either from the centre 
to the exte rnal regi o ns o r limited to the c entral 
iGonzalesI [1991 [Carollo fc Danzigerl [l99l 
1993; Davis et al.l 1993; Fisher & Illingworthl 



region (e.g. 
Carollo et al. 

Il995t iFisher fc Illingworthl 1 1996lh 

The above arguments suggest the occurrence of changes 
in some fundamental properties of the underlying stellar 
population. In order to investigate the stellar population 
and metallicity gradient in NGC 3607, we first analyse the 
spatial distribution of metal- line strengths, and then apply 
a stellar population analysis. 

Strong absorption lines of neutral iron and sodium, as 
well as Mg2 and TiO bands, are present in the spectra 
of NGC 3607 (Fig. [2]). Particularly conspicuous features 
are the Mg2 A5 i 76 , Fei A52 7o, Fei A53 35, Fei A54 o6, Fei A570 9, 
Fel A 5782, Nal A 5895, and TiOn A 6237- These absorption fea- 
tures present a strong dependence with the age and metal- 
licity of the underlying stellar population (e .g lBica fc A lloin 
119861: IRickes et alj[200l IRickes et al.ll2008D . 

Each equivalent width (EW) has been measured three 
times, allowing for the uncertainties in the continuum level 
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Fig. 2. Spatial extractions for the SW (left panel) and NE 
(right panel) directions normalised at A5870A. The dis- 
tance to the galaxy centre is indicated on each spectrum. 
Except for the bottom ones, the spectra have been shifted 
by an arbitrary constant, for clarity purposes 



definition due to noise in the spectra. This procedure al- 
lowed us to estimate the average value and corresponding 
standard deviation from the average for each measurement. 

Subsequently, the observed values of EW were corrected 
for line broadening due to the stellar velocity di spersion. We 
follow ed the same procedure as described in iRickes et al.l 
( 2008|). The spectrum of the G giant star H R 7333, assumed 
to have zero intrinsic velocity, was broadened with a series 
of Gaussian filters with velocity dispersion a varying from 
to 300 km s -1 in steps of lOkrns -1 . For each absorption fea- 
ture considered we calculate an empirical correction index 
C a , so that C<j(EW) = where EW(0) and EW(cr) 

are the observed and broadened equivalent widths, respec- 
tively. We calculate the correction inde x for each extracte d 
spectrum using the a values defined bv lCaon et a l. ( 2000). 

The EWs measured in the spectra of NGC 3607, cor- 
rected for the velocity dispersion, are given in Table [5] 
Except for the Mg2 index, which is measured in magnitude, 
the remaining EWs are given in A. 

The corrected EWs of NGC 3607 (Table HJ) are shown 
as a function of the distance to the galaxy centre in Fig. [3] 
Except for Fe U5782 ; the remaining absorption features 
present significant variations with distance to the centre. In 
particular, the Mg 2 band clearly shows a negative gradient, 



Fig. 3. Radial distribution of the absorption features in 
NGC 3607. The radial scale is given in terms of the effective 
radius, R e — 43.4", with negative values towards the North 
direction. 



i.e. increasing towards the centre of the galaxy. Evidence of 
similar negative gradients are seen in the other Fe I features. 

Similar gradient slopes in the Mg2 index and the 
Fel features suggest the same chemical enhancement pro- 
cess for these elements. Furthermore, as shown in Fig. @] 
the EW of Mg2 correlates both with EW(Fei\ 52 7 ) and 
EW(Fe IA5335), which is consistent with the same enrich- 
ment process acting on Mg and Fe in NGC 3607. 

When the gravitational potential is strong enough to re- 
tain the gas ejected by supernovae events and massive stars 
winds, it eventually migrates to the galaxy's central regions. 
Thus, it should be expected that subsequent generations 
of stars will be more metal rich at the centre than in the 
external parts. Consequently, a negative radial metallicity 
gradient should be established as a mass-dependent param- 
eter. In this scenario, a galaxy formed essentially through a 
monolithic collapse, i.e. with a small contribution of merg- 
ers, is expected to present a well-defined c orrelation be- 
tween the metallicity gradient (iLarsonl 1 19741 : lOgando et alj 
I2OQ50 and total mass (jCarollo et al.lll993f) . 

Th e above issue was investigated by ICarollo et al.l 
(|1993f ). who studied the dependence of the gradient 
<iMg2/<ilog r with galaxy mass in a sample of 42 ellipti- 
cal galaxies. They detected a conspicuous correlation for 
masses lower than ~ 1O 11 M , followed by a saturation 
for higher masses. Different formation mechanisms, such as 
a less dissipative collapse and/or merger of smaller galax- 
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Table 2. Equivalent widths measured in the spectra of NGC3607 



R R/R e Mg 2 index Mg 2 A5176 FeI\ 52 7Q Fel A5335 Fel xswe FeI X570 <j FeI X578 2 Nal\5895 TiOi 



(") 






(mag) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


0.00 


0.00 


0. 


17 ±0.01 


6.10 ±0.22 


3.60 ±0.17 


3.75 ±0.15 


2.27 ±0.04 


1.33 ± 0.02 


1.01 ± 0.05 


5.03 ± 0.21 


6.61 ±0.33 


3.05S 


0.07 


0. 


16 ±0.01 


6.03 ±0.25 


3.58 ±0.35 


3.29 ±0.26 


2.35 ±0.09 


1.34 ± 0.04 


1.06 ±0.01 


5.36 ± 0.22 


6.42 ±0.32 


6.10S 


0.14 


0. 


16 ±0.01 


5.78 ±0.35 


3.46 ± 0.34 


3.13 ±0.52 


2.38 ±0.03 


1.43 ± 0.06 


1.15 ± 0.01 


4.76 ± 0.25 


6.41 ±0.35 


9.15S 


0.21 


0. 


15 ±0.01 


5.49 ±0.38 


3.45 ±0.28 


3.38 ±0.20 


2.08 ±0.11 


1.30 ± 0.03 


1.29 ±0.03 


4.72 ± 0.30 


6.54 ±0.33 


12.20S 


0.28 


0. 


15 ±0.01 


5.56 ±0.38 


3.08 ±0.24 


3.07 ±0.29 


2.02 ±0.15 


1.02 ± 0.10 


1.26 ± 0.08 


4.04 ±0.27 


5.64 ±0.35 


15.25S 


0.35 


0. 


14 ±0.01 


4.97 ±0.26 


2.80 ±0.48 


2.76 ±0.16 


2.35 ±0.06 


1.08 ± 0.08 


1.12 ± 0.02 


4.00 ± 0.24 


5.65 ±0.38 


18.30S 


0.42 


0. 


14 ±0.01 


4.95 ±0.40 


2.80 ±0.56 


2.73 ± 0.30 


2.27 ±0.07 


1.06 ± 0.11 


1.01 ±0.07 


3.99 ± 0.31 


4.85 ± 0.34 


24.40S 


0.56 


0. 


13 ±0.01 


4.89 ±0.53 


3.14 ±0.49 


2.76 ±0.27 


1.86 ±0.14 


1.05 ± 0.11 


1.33 ±0.13 


3.93 ± 0.32 


5.63 ±0.37 


30.50S 


0.70 


0. 


14 ±0.01 


4.92 ± 0.62 


3.21 ±0.63 


2.65 ±0.57 


1.79 ±0.18 


1.23 ±0.12 


0.92 ± 0.02 


3.43 ± 0.32 


5.35 ±0.38 



3.05N 0.07 0.16 ±0.01 6.03 ±0.20 3.64 ± 0.33 3.43 ±0.17 2.28 ± 0.04 1.31 ±0.04 1.00 ±0.02 5.05 ± 0.23 6.38 ± 0.33 

6. ION 0.14 0.16 ±0.01 5.89 ±0.28 3.54 ± 0.37 3.34 ± 0.26 2.26 ±0.15 1.22 ±0.05 1.03 ± 0.03 4.67 ± 0.30 6.36 ± 0.35 

9.15N 0.21 0.15 ±0.01 5.52 ±0.08 3.45 ± 0.24 3.15 ±0.32 2.21 ±0.16 1.30 ±0.04 1.19 ± 0.02 4.96 ± 0.27 6.40 ± 0.31 

12.20N 0.28 0.14 ±0.01 5.39 ±0.28 3.18 ±0.27 3.00 ± 0.24 2.06 ± 0.09 1.07 ±0.07 1.18 ± 0.02 4.56 ± 0.31 5.55 ± 0.33 

15.25N 0.35 0.14 ±0.01 5.09 ±0.16 2.84 ±0.51 2.55 ± 0.24 2.22 ±0.10 1.11 ±0.12 1.03 ± 0.05 4.08 ± 0.25 5.46 ± 0.38 

18.30N 0.42 0.15 ±0.01 5.00 ±0.30 2.81 ± 0.45 2.74 ±0.13 2.21 ±0.10 1.03 ±0.14 0.95 ± 0.02 4.02 ± 0.27 4.78 ± 0.35 

24.40N 0.56 0.14 ±0.01 4.93 ±0.23 2.93 ± 0.50 2.74 ± 0.26 1.93 ±0.13 1.02 ±0.05 1.29 ± 0.05 3.89 ± 0.30 5.65 ± 0.37 

30.50N 0.70 0.13 ±0.01 4.89 ±0.12 3.10 ±0.60 2.59 ± 0.46 1.82 ±0.23 1.19 ±0.10 0.98 ± 0.01 3.55 ± 0.31 5.47 ±0.35 

Table Notes. Note that EW of Mg 2 index is given in mag. Extraction area of all spectra of NGC 3607 is 4.57(") 2 ~ 0.02 kpc 2 . The 
EWs are corrected for velocity dispersion and the errors are the standard deviation of the average of the three different EW 
measurements for each line. 



ies, are invoked to explain the lack of Mgi gradient slope 
to g alaxy mass correlati on for the high-galaxy mass range 
(e.g. ICarollo et alj [19930 . In addition, the correlation be- 
tween e?Mg 2 jd log r with velocity d ispersion (a) - confin ed 
to the low-mass galaxies found by lOgando et al.1 (|2005h - 
suggest that only galaxies below some limiting a have been 
formed by collapse, whereas the massive ones were formed 
predominantly by mergers. 

In order to check how NGC 36 07 fits in the metallicity 
gradient vs. galaxy mass picture of lCarollo et al.l (|1993h . we 
first estim ate the galaxy's dynamic al mass with the approx- 
imation of iPoveda fe Alien! (|1975f ). M tot = 0.9a 2 R e /G = 
3.34 x 1O 8 M0, where G is the gravitational constant, R e = 
28.53" is the effective radius and a — 225 km s -1 is the cen- 
tral stellar velocity dispersion. The dynamical mass com- 
puted for NGC 3607 is M = (1.6±0.3) x 10 11 M which, to- 
gether with the Mg 2 index gradient slope dMg% /d log r = 
—0.02, pu ts this galaxy in the merger-dominated region (see 
Fig. 11 of lCarollo et"aLlfl99l . 

4. Stellar population synthesis 

In this section we investigate the star formation history and 
derive properties of the stellar population of NGC 36 07 by 
means of the stellar population synthesis method of iBical 
(|1988f ) . As popul ation templates we use the synthetic star 
cluster spectra of iBruzual & Charlotl (|2003m 1. 

As a first approach we tested large bases with templates 
with ages between 1 Myr and 13Gyr, but we found that 
the flux contribution of the components younger than a 
few 10 2 Myr is statistically meaningless. Thus, for simplic- 
ity we built a spectral base containing 7 templates. After 
several tests we verified that the base that best reproduces 
the observed spectra, without redundancy, is the one with 

2 Built with a lSafpeterl ljl955h Initial Mass Function 



Table 3. Age and metallicity components 





Z = 0.05 


Z = 0.02 


Z = 0.008 


13Gyr 


Al 


A2 


A3 


5Gyr 


Bl 


B2 


B3 


lGyr 




C2 





the ages 1, 5 and 13Gyr, together with the metallicities 
Z = 0.008, 0.02, and 0.05. Such values are taken as rep- 
resentative of the sub-solar, solar an d above-solar metal- 
licity ranges (e.g. iRickes et al.ll2008lh The distribution of 
templates in the age- metallicity plane is shown in Table [31 
We note that the working base is rather restricted in com- 
pone nts to minimise po tential age/metallicity degeneracies 
(e.g. IRickes et al.l feOOS') due to the relatively small number 
of absorption features available in our spectra of NGC 3607 
(Fig.©. 

We now compare the present stellar pop ulation synthe- 
sis res ults with previous works. Recentelly, lAnnibali et al.l 
(2007j) have studied the underlying stellar population of a 
large sample of early-type galaxies, among them NGC 3607, 
by means of the H/3, Iron and Mgb age and metallicity 
indicators. They report for this galaxy an average stellar 
population age of 3.1 Gyr, while we found an average age of 
8.8 Gyr, based on the the central and external region values, 
a value that is more consistent for an SO galaxy. This dis- 
crepancy could be due to the different methodologies used 
to obtain the age of NGC 3607. Interestingl y, we note that 
this d oes not occur with other galaxies of lAnnibali et all 
(|2007l ) that have been studied by us following the same 
techniques as in the present pape r. For NGC 5044 and 
NGC 6868 (|Rickes et al.ll2004 |2008t) we obtained the aver- 
age ages of 13 and 9.8 Gyr, respec tively, which are very c lose 
to the 14.2 and 9.2 Gyr found by lAnnibali et alj (|2007l h 
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Fig. 4. Correlations among the EWs of selected absorption Fig 5 F i ux _f ract ion contribution at 5870 A of the stellar- 



features measured in NGC 3607. The correlation coefficients 
(CC) of the fits are indicated 



population templates to the spectra of NGC 3607. Negative 
values towards the North direction. 



The stellar population synthesis algorithm provides di- 
rectly the flux fractions (relative to the flux at A5870A) 
that each age and metallicity template contributes to the 
observed spectrum. The sum of the template spectra, ac- 
cording to the individual flux fractions, should be repre- 
sentative of the stellar population contribution to the ob- 
served spec trum. A detailed d escription of the algorithm is 
provided in lRickes et al.|[2004l 

Results of the stellar population synthesis are sum- 
marised in Fig. which shows the 1, 5, and 13 Gyr compo- 
nent contributions to the flux at 5870 A, along most of the 
galaxy's radial extent. The 5 and 13 Gyr curves in Fig. [5] 
correspond to the sum over the three metallicity compo- 
nents. In all cases the solar metallicity is the dominant 
component with a contribution higher than 70%, while the 
Z = 0.05 component contributes with w 10%, and the 
Z = 0.008 with the remaining » 20%. Within uncertainties, 
the 5 and 13 Gyr components dominate the spectra, from 
the centre to the external parts. The error bars in Fig. [5] 
are the standard deviation from the average of all possible 
solutions for each stellar population template. 

As a caveat we note that our observed spectra of 
NGC 3607 cover a relatively short spectral range. Thus, 
they contain a restricted number of indices sensitive to 
age and metallicity that, in principle, could provide a small 
number of constraints for the synthesis, especially in terms 
of metallicity. Nevertheless, Fig [5] indicates that the 13 Gyr 
population dominates in NGC 3607, contributing with ap- 



proximately 50% of the flux at 5870 A. The 5 Gyr popula- 
tion contributes with approximately 30%. 

Another way to assess the efficiency of the stellar popu- 
lation synthesis algorithm is by means of the residuals (dif- 
ference between the observed and synthesised values) , both 
in EWs and continuum points. With a few exceptions, the 
synthesised EWs are within ±lA of the observed values, 
while for the synthesised continuum points (normalised to 
A5870A) the offset lies within ±0.1. We illustrate this in 
Fig. [6] for representative features and continuum points. 
The residuals are almost equally distributed between pos- 
itive and negative values, which indicates that metallicity 
variations in NGC 3607 are essentially contemplated by the 
adopted templates (Table [3]). 

5. Metallicity 

We also investigate metallicity properties of NGC 3607 by 
comparing the measured Lick indices of Mg 2 , Fe IA5270 , and 
Fe IA5335 (Table [2]), with the equivalent one s derived from 
single-aged stellar populat ion (SSP) models flThomas et alJ 
2003, iBuzzoni et al.1 Il992f) that assume a ISalpeterl (|l955h 
initial mass function and different ages. After testing a 
range of SSP ages, we found that the indices measured in 
NGC 3607 are best described by the 10 Gyr models. 

In Fig. [7] we compare our Lick indices measured in 
NGC 3607 with the equivalent ones computed from SSP 
models built for different metallicities and [a/Fe] ratios. 
The central parts of NGC 3607 present a small deficiency of 
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EW A (synt)) of representative EWs (top panel) and contiri- of mctallicity and r a/Fe ] ratios . Filled symbols indicate the 
uum points (bottom). Negative values towards the North centr&1 (higher yalues) and external regions (lower values) 



direction. 



respectively. 



alpha elements with respect to iron (—0.3 [a/Fe] 0.0), 
whereas the opposite trend appears to characterise the 
external parts. We found the average value of [a/Fe] = 
0.13 , which is comp arable to the value 0.24 reported by 
lAnnibali et"aT1 (|2007f ). 

NGC 3607 presents sub-solar metallicity, with the ex- 
ternal parts being more metal-poor than the nucleus. The 
differences observed in [a/Fe] imply different mechanism of 
chemical enrichment for the central and external regions of 
NGC 3607. 



6. Ionised gas 

Emission gas has been detected in a large number of 
elliptical and lenticu l ar galaxies (I Phillips et all Il986t 
iMacchetto et al.l . ll996tlRickes et al.l . l2004l2008[ ).' However, 
the origin of this gas and the ionisation source are not 
yet conclusively established, ft is clear from Fig. [2] that 
NGC 3607 presents conspicuous emission lines throughout 
its body. In this section we investigate the physical condi- 
tions and ionisation source of the gas. 

The properties of the emission gas are derived from line 
fluxes measured in the stellar population subtracted spec- 
tra (Sec H|). The emission line fluxes of Ha, [Nii]a6548j 
[Nii]a6584, [Sii]a6717 and [Sii]a6731 were measured with 
Gaussian functions fitted to the profiles. They are listed 
in Tab. GO 



The spatial distributions of the measured ratios 
[Nil]/ Ha and [SI I}/ Ha are shown in Fig. E While the 
former presents significant changes with distance to the cen- 
tre, the latter is more uniformly distributed. 

Assuming that the emission lines are formed essentially 
by recombination, the number of ionising photons (Q(H)) 

can be calculated as Q(H) = 2^\h^t) (lOsterbrockl 
1989), where as (H°,T) is the total recombination coeffi- 
cient, and UHa [H ,T) is the recombination coefficient for 
Ha. The values of Q(H ) and Lhu are given in Tab. El and 
their spatial distributions are shown in panels (a) and (b) 
of Fig. M 

The central region of NGC 3607 presents [Nil] /Ha = 
3.44 and [Sn]/ Ha > 1. We note that the ratio [Nil]/ Ha is 
higher than 1.0 for all sampled regions (panel (a)). 

Three hypotheses for the nature of the ionising source 
in the external regions of NGC 3607 are tested: the gas is 
ionised by (i) an Hn region, (ii) hot, post-AGB stars, or 
(Hi) an AGN. They are discussed individually below. 



6.1. H II region 

We build an H II region with a mass distributi on that fol- 
lows the initial mass function of Salpcter] (|1955f ) . for stellar 
masses in the range 0.1 — 100 M Q . The model which best 
reproduces the number of ionising photons measured in the 
region at 3.05" is described in Tab. [51 This table provides 
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Table 4. Emission line parameters of NGC 3607 
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R (") 


Ha 


[AT / 6 584] 


[7V// 6 548l 


[Sllmvr] 


\SIhni\ 


[Nll\ 

Ha 


{8U\ 

Ha 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


0.00 


1.34 ± 0.1 


3.46 ± 0.1 


1.15 ± 0.1 


0.66 ± 0.1 


0.69 ± 0.1 


3.44 ± 0.1 


1.01 ± 0.1 


3.05S 


1.29 ± 0.1 


3.20 ± 0.1 


1.06 ± 0.1 


0.58 ± 0.2 


0.72 ± 0.2 


3.30 ± 0.1 


1.01 ± 0.1 


6.10S 


0.93 ± 0.1 


2.20 ±0.1 


0.73 ±0.1 


0.47 ±0.2 


0.50 ±0.2 


3.15 ±0.1 


1.04 ± 0.1 


9.15S 


0.60 ± 0.2 


1.43 ± 0.2 


0.47 ± 0.2 


0.32 ± 0.2 


0.28 ± 0.2 


3.17 ± 0.1 


1.00 ± 0.1 


12.20S 


0.56 ± 0.2 


1.10 ± 0.2 


0.36 ± 0.2 


0.14 ± 0.1 


0.32 ± 0.1 


2.61 ± 0.1 


0.82 ± 0.1 


15.25S 


0.60 ± 0.2 


1.02 ±0.2 


0.34 ±0.2 


0.22 ±0.2 


0.35 ±0.2 


2.27 ±0.1 


0.95 ± 0.1 


18.30S 


0.84 ±0.2 


1.25 ±0.2 


0.41 ±0.2 


0.19 ±0.2 


0.36 ±0.2 


1.98 ±0.1 


0.65 ± 0.1 


24.40S 


0.80 ± 0.2 


1.07 ±0.3 


0.35 ±0.2 


0.28 ±0.2 


0.30 ±0.2 


1.78 ±0.1 


0.73 ± 0.1 


30.50S 


0.85 ± 0.3 


0.89 ± 0.3 


0.29 ± 0.3 


0.26 ± 0.2 


0.31 ± 0.2 


1.39 ± 0.1 


0.67 ± 0.1 


3.05N 


0.97 ±0.1 


2.43 ±0.1 


0.81 ±0.1 


0.42 ±0.2 


0.44 ±0.2 


3.34 ±0.1 


0.89 ±0.1 


6.10N 


0.71 ± 0.1 


1.87 ±0.1 


0.62 ±0.1 


0.41 ± 0.2 


0.47 ±0.2 


3.51 ±0.1 


1.24 ± 0.1 


9.15N 


0.66 ± 0.2 


1.43 ±0.2 


0.47 ±0.2 


0.31 ± 0.2 


0.31 ±0.2 


2.88 ±0.1 


0.94 ± 0.1 


12.20N 


0.65 ± 0.2 


1.16 ±0.2 


0.38 ±0.2 


0.26 ± 0.2 


0.26 ±0.2 


2.37 ±0.1 


0.80 ±0.1 


15.25N 


0.80 ± 0.2 


1.07 ±0.3 


0.35 ±0.2 


0.30 ± 0.2 


0.30 ±0.2 


1.78 ±0.1 


0.75 ±0.1 


18.30N 


0.90 ± 0.3 


1.27 ±0.3 


0.42 ±0.3 


0.22 ±0.3 


0.22 ±0.3 


1.88 ±0.1 


0.49 ± 0.1 


24.40N 


0.51 ±0.4 


1.03 ±0.4 


0.34 ±0.4 


0.42 ±0.3 


0.42 ±0.3 


2.69 ±0.4 


1.65 ± 0.2 
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Table 5. Ha luminosity and number of ionising photons 
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Fig. 8. Spatial variation of the ratios [JV/ ^ 6584 (a) and 
^ SII Ha™ 0-0 ■ Negative values towards the North direction. 



the distribution of stars for a representative set of spectral 
type and mass; in the last column we give the number of ion - 
ising photons (A^y a ) produced by each star (|Dottorilll98ll) . 
The total number of ionising photons can be computed from 



R 


F Ha 


LHa 


Q(H) 


(") 


(10" 15 erg cm- 2 s _1 ) 


(10 37 ergs _1 ) 


(10 49 s" 1 ) 


o.oos 


1.34 ±0.10 


2.61 ±0.03 


13.7 ±0.02 


3.05S 


1.29 ± 0.10 


2.51 ± 0.02 


13.2 ± 0.02 


6.10S 


0.93 ±0.01 


1.81 ±0.03 


9.54 ±0.01 


9.15S 


0.60 ± 0.02 


1.17 ±0.02 


6.15 ±0.02 


12.20S 


0.56 ± 0.02 


1.09 ± 0.02 


5.74 ± 0.03 


15.25S 


0.60 ± 0.02 


1.17 ±0.03 


6.15 ±0.03 


18.30S 


0.84 ±0.02 


1.63 ±0.02 


8.62 ±0.03 


24.40S 


0.80 ± 0.02 


1.56 ±0.02 


8.21 ±0.01 


30.50S 


0.85 ± 0.03 


1.65 ±0.02 


8.72 ± 0.02 


3.05N 


0.97 ±0.01 


1.89 ±0.03 


9.95 ± 0.02 


6.10N 


0.71 ±0.01 


1.38 ±0.03 


7.28 ±0.01 


9.15N 


0.66 ± 0.02 


1.28 ± 0.03 


6.77 ± 0.01 


12.20N 


0.65 ± 0.02 


1.26 ± 0.03 


6.67 ±0.02 


15.25N 


0.80 ± 0.02 


1.56 ± 0.03 


8.21 ± 0.02 


18.30N 


0.90 ± 0.03 


1.75 ±0.03 


9.23 ± 0.02 


24.40N 


0.51 ±0.04 


0.99 ± 0.03 


5.23 ±0.01 



Table Notes. Ha 
computed for 



N(Ly) = N Ly (m)(j)(m)dm. The relation between the 
number of ionising photons with stellar mass can be approx- 

/ \5.9 

imated as N Ly (m) w 3.9 x 10 39 f ^-j s' 1 (|Rickes et all 

2008) . Given the above, the adopted model H II region pro- 
duces ~ 2.6 x 10 50 photons s" 1 , which is comparable to the 
observed value. However, with more than 20 O stars, such 
an H II region would be very young, massive, and luminous 
enough to have already been detected in previous works. 
Besides, the stellar population synthesis (Sect.0]) does not 
indicate significant, recent star formation in the sampled 
regions of NGC 3607. 



luminosity and 
the extractions e 



number of ionising photons 
iven in col. 1. 
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Fig. 9. Spatial distribution of the Ha luminosity (a) flux 
(c) and number of ionising photons (b). Negative values 
towards the North direction. 

Table 6. Hn region model parameters 



Spectral type 


M 


Number 


N Lya 




(M ) 


(stars) 


(10 48 photons" 1 ) 


05 


30 


20 


2.00 




22 


41 


0.25 


B0 


17 


75 


0.082 




15 


101 


0.043 




13 


142 


0.015 




10 


264 


0.0021 



Table Notes. Col. 2: representative stellar mass; Col. 3: number 
of stars; Col. 3: number of ionising photons. 



6.2. Post- AG B stars 

The number of ionising ph otons emitted by a typical Post- 
AGB star is - lO^s" 1 (jBinettd . fl99l . Thus, about 10 4 
such stars would be necessary to produce the amount of 
photons measured in each extracted region. Consequently, 
for an extraction area of « 0.02 kpc 2 (Sect. [2]), the pro- 
jected surface-density of post-AGB stars would b e about 
5 x 1 5 kpc~ 2 , a number prohibitively high (e.g. iBinettd 
119941) . Besides, the spatial profiles of [Nil] and Ha ex- 
tracted along the slit (Fig. [8j, show a ratio N[n]/Ha > 
2 for most regions, especially towards the central parts. 
However, the relatively high emission-line ratios measured 
in NGC 3607, especially in the central region, are not con- 
sistent with ionisation models based on hot old stars in el- 



liptica l galaxies (e.g. lBinettdll994f ). The models of lBinette] 
(|1994f ) do not reproduce as well the observed large ratios 
[Nil] /Ha and [S ii] 6 y 31 /Ha. Instead, high values for these 
emission-line ratios are typical of AGNs. 

On the other hand, the lower ratios of [N 1 1] A6584 /H a 
and [SI I] 6731 /H a (Fig. [8]) together with the lower values 
of Q(H) (Fig. measured for large radii may, in part, be 
accounted for by hot old stars. 

6.3. Active nucleus 

The arguments that follow from the above discussion sug- 
gest that NGC 3607 hosts an AGN. We test this hypothe- 
sis with the emission-line diagnostic-diagram involving the 
ratios [NII]\ 65S 4/Ha and [S'//]a6717,3i/ Ha. For compar- 
ison we take the equivalent values measured in a sam- 
ple of nearby galaxies that host a low-luminosit y (Lh q < 
2 x lO^ergs" 1 ) active nucleus (jHo et al.|[T997t) . To min- 
imise the possibility of uncer tain classi f icatio n, we restrict 
the analysis to the galaxies in lHo et"aH (|1997f ) sample with 
unambiguous classification as LINER in spiral, LINER in 
elliptical, Seyfert or S tarburst. The restri cted sample con- 
tains 128 galaxies (see lRickes et al.lfeOOSD . 

The results are shown in Fig. [TfJl where it can be seen 
that the values measured in NGC 3607 occur in the locus 
occupied predominantly by LINERs and Seyferts. However, 
the lack of a broad component in Ha (Fig. \2§ strongly 
favours the LINER interpretation. Together with the stel- 
lar population synthesis, this result suggests that the main 
source of gas ionisation in NGC 3607 is non-thermal, pro- 
duced by a low-luminosity active nucleus. The lower ratios 
of [NII]\ 65S 4/Ha and [SI 7] 6731 /Ha towards the external 
parts can be accounted for by both a LINER and hot stars. 

Finally, in Fig.[TTJwe test whether the gradients in emis- 
sion and absorption-line features are somehow related. We 
work with the EW of Mg 2 A5176 (Table C2) and the Ha 
flux (Table |4j distributions, taking the average values of 
the North and South extractions (the 2 outermost points 
have been averaged out to minimise error bars). Within 
the uncertainties, both distributions appear to correlate for 
R Si 12", which can be accounted for the steep and mono- 
tonic Ha flux increase towards the central region (Fig. [5]) 
together with the overall negative Mg2 gradient (Fig. [3]). 
While the EW of Mg 2 keeps dropping for R ^ 12", the 
Ha flux oscillates abruptly, probably because of the rela- 
tive contribution of the hot, old stars to the gas ionisation, 
superimposed on the severely diluted LINER flux at such 
large radii. This combination of effects is reflected on the 
mild anti-correlation of both properties (lower portion of 
Fig. ITT]) . It is also consistent with the low ratios of [N 11] /H a 
and [Sn]/Ha (Fig. [5]), and the low Q(H) (Fig.©, measured 
for large radii, all accounted for by hot, old stars. 

Similar results would be obtained with other absorption 
features, since they also present similar gradients (Fig. [3]). 

7. Discussion and concluding remarks 

Individual galaxies can provide important pieces of evi- 
dence that may ultimately lead to a better understand- 
ing of the formation and evolution processes of the early- 
type galaxies as a class. Following the methods previously 
employed in the analysis of other ea rly- type galaxies (e.g. 
iRickes et al1l2004HRickes et al.ll2004lh in the present paper 
we focus on the lenticular galaxy NGC 3607. 
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Fig. 10. Emission-li ne diagnostic di agrams with a subsam- 
ple of the galaxies in lHo et al. ( 1997) plotted as comparison. 
The locus of NGC 3607 is consistent with that of LINERs. 



We investigate the stellar population, metallicity dis- 
tribution and ionised gas of NGC 3607 by means of long- 
slit spectroscopy (covering the wavelength range AA5100 — 
6800A) and a previously tested stellar population synthe- 
sis algorithm. The stellar population synthesis shows dif- 
ferent star formation episodes with ages distributed within 
1-13 Gyr. 

The radial distribution of the equivalent widths mea- 
sured for the available metallic absorption features in this 
galaxy consistently present a negative gradient, which in- 
dicates an overabundance of Fe, Mg, Na and TiO in the 
central parts with respect to the external regions. Also, the 
central region of NGC 3607 presents a deficiency of alpha 
elements and higher metallicity than the external parts, 
which suggests different star-formation histories in both 
regions. With respect to the emission gas, the line ratios 
[NII)\ 658 ,4/Ha and [SII]x67n,3i/ Ha consistently indicate 
the presence of a LINER at the centre of NGC 3607. 

We found that the Mg^ index correlates both with 
FelA527o an d FelA5335 ■ which suggests that Mg and Fe prob- 
ably underwent the same enrichment process in NGC 3607. 
Also, the Ha flux distribution correlates with the radial 
profile of metallic features, in the sense that they both de- 
crease from the galaxy centre to R sa 12". The correlation 
breaks for larger radii, especially because, while the metallic 
feature profiles keep dropping, gas ionisation due to hot, old 
stars adds to the diluted LINER flux. Considering the re- 
lation between galaxy dynamical mass an d the Mg2 gradi- 
ent slope (Figure 11 of lCarollo et al.|[l993l ). the parameters 
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Fig. 11. Relation between EW(Mg 2 A5176) and the Ha 
flux as a function of distance to the galaxy centre. 

computed for NGC 3607 in the present work suggest that 
this galaxy may have underwent important merger events 
in the past. 
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